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The unconstrained classical system equivalent to spatially homogeneous SU (2) Yang- 
Mills theory with theta angle is obtained and canonically quantized. The Schrodinger 
eigenvalue problem is solved approximately for the low lying states using variational cal- 
culation. The properties of the groundstate are discussed, in particular its electric and 
magnetic properties, and the value of the "gluon condensate" is calculated. Furthermore 
it is shown that the energy spectrum of SU{2) Yang-Mills quantum mechanics is inde- 
pendent of the theta angle. Explicit evaluation of the Witten formula for the topological 
. . . susceptibility gives strong support for the consistency of the variational results obtained. 
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. I. INTRODUCTION 

> 

o 

f — . ' One of the central problems in the investigation of non-Abelian gauge theories is a gauge invariant 
, description of the vacuum and the low lying exited states. In the standard approach to the quantization 
of gauge theories the physical states have to satisfy not only the Schrodinger equation but additionally be 
unnihilated by the Gauss law operator to implement gauge invariance at the quantum level ||l| . However 
it is well known that there exist states which satisfy the Gauss law but are not invariant under the so- 
called homotopically nontrivial gauge transformations, leading to the appearence of the theta angle 
G |§| • A well-elaborated semiclassical approach to the theta structure of the groundstate has been given in 
I ' the "instanton picture" , where the theta angle is interpreted 0] in analogy to the Bloch momentum in 
Qh, Solid State Physics. The instantons, which are selfdual solutions of the Euclidean classical equations of 
^ ' motion with finite action, correspond to semiclassical quantum mechanical tunneling paths in Minkowski 
space between the infinite sequence of degenerate zero-energy Yang-Mills vacua of different homotopy 
classes of the gauge potential. The semiclassical instanton picture of the theta vacuum however is of 
course reliable only for weak coupling. For a complete investigation of the theta stucture of the vacuum 
of Yang-Mills quantum theory a rigorous treatment at strong coupling is necessary. The effect of the 
■ theta angle for arbitrary coupling constant can be taken into account by adding the Pontryagin density 
to the Yang-Mills Lagrangian [Q . Although the extra theta dependent CP-violating term is only a total 
divergence and therefore has no meaning classically, it can have a physical meaning at the quantum level 
as is still under lively discussion [Q- 

As a first step towards a full investigation of Yang-Mills theory in the strong coupling limit the toy 
model of SU{2) Yang-Mills mechanics of spatially homogeneous fields has been considered on the classical 
- 1 11 1 as well as on the quantum level [|l2|- |^l[. In the present paper we will analyse the model 



of SU(2) Yang-Mills mechanics of spatially homogeneous fields for arbitrary theta angle. In order to 
obtain the equivalent unconstrained classical system in terms of gauge invariant variables only - 
p5[ , we apply the method of Hamiltonian reduction ( | |35| and references therein) in the framework of 
the Dirac constraint formalism - ||3^ . As in our recent work jl|] the elimination of the pure gauge 
degrees of freedom is achieved by using the polar representation for the gauge potential, which trivializes 
the Abelianization of the Gauss law constraints, and finally projecting onto the constraint shell. The 
obtained unconstrained system then describes the dynamics of a symmetrical second rank tensor under 
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spatial rotations. The main-axis-transformation of this symmetric tensor allows us to separate the gauge 
invariant variables into scalars under ordinary space rotations and into "rotational" degrees of freedom. In 
this final form the physical Hamiltonian and the topological operator can be quantized without operator 
ordering ambiguities. We study the residual symmetries of the resulting unconstrained quantum theory 
with arbitrary theta angle and reduced the eigenvalue problem of the Hamiltonian to the corresponding 
problem with zero theta angle. Using the variational approach we calculate the low energy spectrum with 
rather high accuracy. In particular we find the energy eigenvalue and the magnetic and electric properties 
of the groundstate, as well the corresponding value of the "gluon condensate". The groundstate energy 
is found to be independent of the theta angle by construction of the explicit transformation relating 
the Hamiltonians with different theta parameter. This is confirmed by an explicit calculation of the 
Witten formula for the topological susceptibility using our variational results for the groundstate and the 
low lying excitations give strong support for the independence of the groundstate energy of theta thus 
indicating the consistence of our results. 

Our paper is organized as follows. In Section II the Hamiltonian reduction of SU{2) Yang-Mills 
mechanics for arbitrary theta angle is carried out and the corresponding unconstrained system put into 
a form where the rotational and the scalar degrees of freedom are maximally separated. In Section HI 
the obtained unconstrained classical Hamiltonian is quantized, its residual symmetries, the necessary 
boundary conditions for the wave functions and the relevance of the theta angle on the quantum level 
discussed. In Section IV the eigenvalue problem of the unconstrained Hamiltonian with vanishing theta 
angle is solved approximately in the low energy region using the variational approach. In Section V 
the Witten formula for the topological susceptibility is evaluated using the obtained variational results. 
Section VI finally gives our conclusions. Appendices A to C state several results and additional discussions 
relevant to the main text. 



II. UNCONSTRAINED CLASSICAL SU{2) YANG-MILLS MECHANICS WITH THETA 

ANGLE 



A. Hamiltonian formulation 



It is well known 0] that the theta angle can be included already at the level of the classical action 



S[A] ■■=-{! fc.F'^'^'^ - ^F;:,F^'^-^ ) , (2.1) 



with the SU{2) Yang-Mills field strengths F^^^ df,A'i;-d,,A1^+ge'''"'A''^A''^ , (a 1, 2, 3), the dual Fa^^ 
l/2e u,,ya-oF°''^P and as = /An. For the special case of spatially homogeneous fields the Lagrangian in 
(2.1) reduces to [] 

i = 2 ^ g^ahcAhQAcij - -Bli - yAai - geabcAboAdj Bai , (2.2) 

with the magnetic field Bai — {^/2)g€abc^ijkAbjAck- After the s upp osition of spatial homogeneity of the 
fields the SU{2) gauge invariance of the Yang-Mills action Eq. (2.1) reduces to the symmetry under the 
SO{3) local transformations 

Aaoit) ^ A^oit) = OiLo{t))abAboit) - ^eabc [o{uj{t))d{io{t)) 

Aa^it) ^ A^^^it) - OiL0{t))abAUt) (2.3) 

and as a result the Lagrangian (^.2|) is degenerate. From the calculation of the canonical momenta 



^ Everywhere in the paper we put the spatial volume V = 1. As result the coupling constant g becomes 
dimensionful with g^^^ having the dimension of energy. The volume dependence can be restored in the final 
results by replacing g^ by g^ /V. 
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Pa := dL/didoAao) = , n,, := dL/d{doAai) = A„ ~ geabcAmA,^ - , (2.4) 

ZTT 

one finds that the phase space spanned by the variables {AaQ^Pa) and [Aai^I^ai) is restricted by the 
primary constraints Pa{x) = . The evolution of the system is governed by the total Hamiltonian |Q 
with three arbitrary functions \a{x) 



Ht ■■= t;K, + o 1 + B^A) + OQiU, A) ~ gAaoeabcAuTl,, + Xaix)Paix) , (2.5) 




where the topological charge has been introduced 

(3:--^n„B,, . (2.6) 

ZTT 

Apart from the primary constraints P" = the phase space is restricted also by the non-Abelian Gauss 
law, the secondary constraints 

$a := geabcAcUM = , = ge^fe^fc , (2.7) 

which follow from the maintenance of the primary constraints in time. 

To overcome the problems of the existence of these constraints and the nonuniquc character of the 
dynamics governed by the total Hamiltonian ( 2^ ) we will follow the method of Hamiltonian reduction 
to constuct the unconstrained system with uniquely predictable dynamics. As in the recent paper p| 
we shall use a special set of of coordinates which is very suitable for the implementation of Gauss law 
constraints and the derivation of the physically relevant theory equivalent to the initial degenerate theory. 
This will be the subject of the following Subsection. 



B. Canonical transformation to adapted coordinates and projection to Gauss law constraint 



The local symmetry transformation (2.2) of the gauge potentials Aai promts us with the set of coor- 
dinates in terms of which the separation of the gauge degrees of freedom occurs. As in we use the 
polar decomposition for arbitrary 3x3 quadratic matrices [pa| 



Aai iX, S) = Oak (X) Skt , 



(2.8) 



with the orthogonal matrix 0{x), parametrized by the three angles Xi ^^d the positive definite 3 x 
3 symmetric matrix S. The representation (2.8) can be regarded as transformation from the gauge 
potentials Aai to a the set of coordinates Xi and 5*^^. The corresponding canonical conjugate momenta 
{pxi, Pik) can be obtained using the generating function 



F (H; X,S)=Y^ IVa.Aa, (x, S) = tr {U^O{x)S) 



(2.9) 



as 



Pxi " "n ~ / , ^ai ~ Ssi — tr 

dXj ^ 



c>Xj 



P,.k = ^ = l{OU^+IlO^: 



ik 



ik 



(2.10) 



(2.11) 



A straightforward calculation yields the following expressions for the field strengths Hai in terms of 
the new canonical variables 



Hai = Oakix) 



Pki+£kH{s ^{x)Px).j - '^ri 



(2.12) 



with 



3 



^ij (x) 2 



dO^ ix) 
9Xj 



0{x) 



and 



Sik '■= Sik — SiktrS . 



(2.13) 



(2.14) 



Using the representations ( 2^ ) and ( |2.12| ) one can eas ily convince oneself that the variables S and P 
make no contribution to the Gauss law constraints (2.7) 



(2.15) 



Hence, assuming the invertibility of the matrix Q, the non-Abelian Gauss law constraints are equivalent 
to the set of Abelian constraints 



(2.16) 



After having rewritt en th e model in terms of adapted canonical pairs and after Abelianization of the 
Gauss law constraints (2.7) the construction of the unconstrained Hamiltonian system can be obtained 
as follows. The physical unconstrained Hamiltonian, defined as 



Hg{S,P) Ht{S,P) 



takes the form 



He - -tr(£2) + ^(1 + 002) [tr^iSr - tr(5)4] + 0QiS, P) . 



where the "physical" electric field strengths £ai are 



and the topological charge 



Ha. =: Oakiq) £m{S,P) , 



Q(5,P) = --^tr(P5) 
zvr 



(2.17) 



(2.18) 



(2.19) 



Using the representation ( 2.12| ) for the electric field one can express the Sai in terms of the physical 
variables P and S 

1 



ek,{S,P) = P,k + ^^{sMs)^^ 



(2.20) 

where A4 denotes the spin part of the angular angular momentum tensor of the initial gauge field 

M„„:=(5F-P5),„„ . (2.21) 



Using ( 2.20 ) the unconstrained Yang-Mills Hamiltonian reads 

He {S,P) = hriPf + ^^tr {sMsf + ^ (l + ^O'^ [tr^Sf tr(5)^] + eQiS,P) . (2.22) 



C. Unconstrained Hamiltonian in terms of rotational and scalar degrees of freedom 



In order to achieve a more transparent form for the reduced Yang-Mills system ( p.22| ) it is convenient 
to decompose the positive definite symmetric matrix S as 

S = R^{a,f3,-/) D{xi,X2,X3) R{a,f3n) , (2-23) 
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with the SO{3) matrix R parametrized by the three Euler angles {a, (3, 7), and the diagonal matrix 



D := diag {xi,X2,xz) 



(2.24) 



Using the Xi and the Euler angles (a,/3, 7) and the corresponding canonical momenta pi and pa,Pf3,p-y 
as the new set of canonical variables on the unconstrained phase space we get the following physical 
Hamiltonian 



He ixi,pi]^i) 



2 ^ 

cyclic 



-it 



2 



[x] 



5^(1 



QQ{p,x) 



In ( 2.25 ) all rotational variables are are combined into the quantities 

sin 7 



sin/3 
cos 7 



sin/? 



Pa + cos ^ pp — sin 7 cot /? , 
Pa + sin 7 p/3 + cos 7 cot /3 , 



P7 



representing the 50(3) invariant Killing vectors with the Poisson brackets algebra 



(2.25) 

(2.26) 

(2.27) 
(2.28) 

(2.29) 



The topological charge Q is independent of the rotational degrees of freedom and depends on the diagonal 
canonical pairs in the particularly simple cyclic form 



= -97r (X1^2P3 + X2X3P1 + X3X1P2) 



(2.30) 



This completes our reduction of the spatially homogeneous constrained Yang-Mills system with theta 
angle to the equivalent unconstrained system describing the dynamics of the physical degrees of freedom. 

If we would restrict our consideration only to the classical level, the above generalization to arbitrary 
theta angle would be unnecessary, because the theta dependence enters the initial Lagrangian in the 
form of a total time derivative and thus the value of the theta angle has no influence on the classical 
equations of motion. In the Hamiltonian formulation one can easily verify that the theta dependence can 
be removed from the Hamiltonian Hg by the canonical transformation to the new variables 



Pt 



g^XjXk, J, fc cyclic. 



(2.31) 



However, the transition to the quantum level requires a more care ful treatment of the problem. It 
is necessary to clarify whether the operator corresponding to (2.31) acting on the quantum states is 
unitary. 

In subsequent Sections we shall consider the quantum treatment of the obtained classical system and 
shall discuss the theta dependence of the vacuum in this model. 



III. QUANTIZATION, SYMMETRIES AND BOUNDARY CONDITIONS 



The Hamilton operator corresponding to ( ^.25 ) is obtained in the Schrodinger configuration represen- 
tation by the conventional representation for the canonical momenta p^. = —id/dx^ 



He 



1 ^ 



cyclic _ 



dx. 



a_2 + ii 



(-1 



2\2 



4) 



47r2 



2 2 

XjXf^ 



(3.1) 



with the topological charge operator 
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^? = *5gE^«^^a^' (3-2) 

cyclic 

and the intrinsic angular momenta ^ obeying the commutation relations 

= -i^ijkik ■ (3.3) 

The transition to the quantum system in this adapted basis is free from operator ordering ambiguities. 
As already mentioned in the last section the parameter theta is unphysical on the classical level, since 



it can be removed from Hamiltonian Hg by the canonical transformation (2.31). One can easily convince 



oneself that the quantum Hamiltonians Hg and Hg^Q can be related to each other via the transformation 

Hg = Ui9)Hg=oU-\e) , (3.4) 

with 

U{9) = exp[ig-^6lxia;2a;3] . (3.5) 

The question is whether this operator is unitary in the domain of definition of the Hamiltonians Hg and 
Hg^Q, which is determined by their respective symmetries and the boundary conditions to be imposed 
on the corresponding wave functions. 



A. Boundary conditions 



Due to the positivity of the coordinates Xi in the polar decomposition ( p.SD the configuration space 
is after the elimination of the pure gauge degrees of freedom. Thus the implementation of the 
canonical rules of quantization to the unconstrained classical system requires the specification of the 
boundary conditions both at positive infinity and on the three boundary planes Xi = , i = 1, 2, 3. The 
requirement of Hermiticity of the Hamiltonian Hg (3.1) leads to the condition 



(■^*gdk<i>e - dk%^e + 2ig^9x^Xj%<^>( 
\ Zn 



— , hj,k cyclic 



(3.6) 



Using the relation Vl/g = U{6)'iig^o with U (9) given in (3.5), this reduces to the corresponding requirement 
for the Hermiticity of i?e=o 



(*e=o5fc*e=o - 9fc*;^o*e=o) = , ij,k cyclic. 

Xk=0 



(3.7) 



It is satisfied for (k arbitrary c-number) 



(afc*e=o + K*e=o) 



Xk=0 



A: = 1,2,3 , 



(3.8) 



which includes the two limiting cases of vanishing wave function (k oo) or vanishing derivative of the 
wave function (k — 0) at the boundary. The requirement of the Hermiticity of the momentum operators 



in the Schrodinger configuration representation pi 
the boundary conditions 



-id/dxi on requires the wave function to obey 



*e=o =0, i= 1,2,3 



0, i = l,2,3 



(3.9) 
(3.10) 



In particular, they also imply the Hermiticity and the existence of a real eigenspectrum of the topological 
charge operator Q. Its eigenstates, however, given explicitly in Appendix A, do not satisfy the boundary 
conditions (3.9) and (B.IO), similar to the eigenstates of the momentum operator ~id/dxi. Furthermore, 
it is interesting to note that the characteristics of the Q operator coincide with the Euclidean self(anti- 
)dual zero-energy solutions of the classical equations of motion. They are the analogs of the instanton 
solutions, but do not correspond to quantum tunnelling between different vacua (see Appendix A). 
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B. Symmetries of the Hamiltonians He and Ho^q 



As a relict of the rotational invariance of the initial gauge field theory the Hamiltonian (3^) possesses 
the symmetry 

[i?,Jfc]=0, (3.11) 

where Ji — Rij^j are the spin part of the generators of the angular momentum of Yang-Mills fields 
satisfying the so(3) algebra 

[Ji,Jj] = i^ijkJk , (3-12) 

and commuting with the intrinsic angular momenta, — . Hence the eigenstates can be classified 

according to the quantum numbers J and M as the eigenvalues of the spin = Jf + J2 + and J3 . 
The Hilbert spaces of states with different spin J are each invariant subspaces under the action of all 
generators Ji and can therefore be considered as separate eigenvalue problems. 

Apart from this continous rotational symmetry the Hamiltonians Hg and Hg^Q possess the following 
discrete symmetries. Both Hg^Q and Q are invariant under arbitrary permutations of any two of the 
variables aijXi — XjOij, aijPi — Pjcrij 

[i/e=o,Ty] -0 , [Q,ay]=0. (3.13) 

However, under time reflections Txi = x^T, Tpi — —piT, as well as under parity reflections Vxi = 
—XiV, Vpi = —piV, Hg^o commutes with T and V, 

[i?e=o,T]=0, [Hg=o,V] = 0, (3.14) 

but Q anticommutes with T and T', 

QT = -TQ, QV = -VQ. (3.15) 

Hence for the Hg^o Schrodinger eigenvalue problem we can restrict to the Hilbert space of real and 
parity odd wave functions which automatically satisfy the boundary conditions (|3.9[). Observe that the 



transformation (3^) leads out of the corresponding Hilbert space and is therefore not unitary. 



C. Independence of the energy spectrum of the theta angle 



Due to the relation ( p.4[ ) between the Hamiltonians Hg , and Hg^Q and the corresponding compatibility 
of the boundary conditions discussed above the energy spectrum should be independent of the theta angle. 
In particular the topological susceptibility of the vacuum should vanish. Using the Witten formula pOj , 
|(41| , the topological susceptibility can be represented as the sum of a propagator term involving the 
transition matrix elements of the topological operator Q and a contact term prportional to the vacuum 
expectation value of the square of the magnetic field. Independence of the groundstate energy of the 
theta angle and hence vanishing topological susceptibility should therefore imply 

d'Eoie) _ ^^|(0|Q|n)| 



_25:»|J! + (0|(fi)^B^|0)=0, (3.16) 



where \n > are eigenstates of the Hamiltonian Hg^o with energy eigenvalues En- As we shall see below 
that our calculation of the low energy part of the spectrum of Hg^Q using the variational technique is in 



full accordance with (3.16). 



IV. SCHRODINGER EIGENVALUE PROBLEM FOR VANISHING THETA 

A. Low energy spin-0 spectrum from variational calculation 

The Hilbert space of states with zero spin f = is an invariant subspace under the action of all 
generators Ji and one can consider the eigenvalue problem separately from states characterized by higher 
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spin value. Thus in the sector of zero spin = ^ = the Schrodinger eigenvalue problem (^]|) reduces 
to 



1 ^ 
- F 

2 ^ 

cyclic 



dx. 



2 2 2 
2 ' 9 Xj^k 



(4.1) 



We shall use the boundary conditions (3.!;) and ( 3.1Cl| ). Already a long time ago it has been proven by F. 
Rellich that Hamiltonians of the type (4.1) have a discrete spectrum due to quantum fluctuations, 
although the classical problem allows for scattering trajectories (see discussion in |12 ). Related and 
simplified versions of the eigenvalue problem ( |4.1| ) have been studied extensively by many authors using 
different methods ||l^- |2^]. In particular, in |14||- |T5| the eigenstates and eigenvalues have been found 
in the semiclassical approximation for the special two dimensional case = 0. ^ 

To obtain the approximate low energy spectrum of the Hamiltonian in the spin-0 sector we will use 
the Rayleigh-Ritz variational method [[46[ based on the minimization of the energy functional 



(^|go|^) 



(4.2) 



The key moment in all variational calculations is the choice of the trial functions. Guided by the harmonic 
oscillator form of the valleys of the potential in (|4.l | ) cl ose to the bottom a simple first choice for a trial 
function compatible with the boundary conditions ( p.9[ ) and ( 3.1C| ) is to use [|l| the lowest state of three 
harmonic quantum oscillators on the positive half line 



000 



i=l 



1/4 



LUiXiC 



(4.3) 



The stationarity conditions for the energy functional of this state, 

3 



E 

cyclic 



9 2 1 



;5 



8 UJjLOk 



lead to the isotropic optimal choice 

uj ■= uji = UJ2 = = 3^/^.g^/^ . 
As a first upper bound for the groundstate energy of the Hamiltonian we therefore find 

E, < £[^ooo] = ?3i/3g2/3 = 4.8676 g^/s. 



(4.4) 



(4.5) 



The upper bound ( |4.5| ) is in agreement with the lower bound of the energy functional for separable 
functions 



£[^sep] > 4.5962 5^/3 ^ 



(4.6) 



derived in Appendix B. 



^ It is interesting that for the three dimensional case one can write the potential term in ( [l.l| ) in the form 
V — "Y^^^iidiW)^ with the "superpotential" W{xi,X2, xg) — X1X2X3. Note that in the simplified two dimensional 
case there is no such superpotential. The two-dimensional superpotential W^-^^ = xy corresponds to the two- 
dimensional harmonic oscillator l/*-^^ — x^ + y^. From the form of the superpotential it follows that the wave 
function 'I'o = exp[— pW] solves the Schrodinger eigenvalue problem with energy eigenvalue E = 0. It is the 
unconstrained, strong coupling form of the well-known exact but nonnormalizable zero-energy solution |^ of the 
Schrodinger equation of Yang-Mills field theory . Obviously it is also not satisfying the boundary conditions ( [3.9| ) , 
(3.10) and has to be disregarded as a false groundstate. 
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In order to improve t he u pper bound for the groundstate energy of the Hamiltonian Hq we extend the 
space of trial functions (O) and consider the Fock space of the orthonormal set of products 

3 

1=1 

of the odd eigenfunctions of the harmonic osciUator 



with the frequency fixed by (4.4). 

Furthermore the variational procedure becomes much more effective, if the space of trial functions is 
deco mposed into the irreducible representations of the residual discrete symmetries of the Hamiltonian 
(4.1). As has been discussed in Section III B, it is invariant under arbitrary permutations of any two of 
the variables aijXi — Xjaij, (TijPi — pjCTij and under time reflections Txi = XiT, Tpi — —piT, 

We shall represent these by the permutation operator ai2, the cyclic permutation operator (T123 and the 
time reflection operator T, whose action on the states is 

a'i23^'(a;i,a:2,a;3) = ^'(2:2, 0:3, xi) , 
ai2^{xi,X2,X3) = ^'(x2,a:i,a;3) , 

T'^{xi,X2,X3) ^^*{xi,X2,X3) , 

and decompose the Fock space spanned by the functions ( [4.7D into the irreducible representations of the 
permutation group and time reflection T. For given (rii, n2, n^) we define 

"f^nlt ■■= *n„n , (4.8) 

if all three indices are equal (type I), the three states (m = —1,0, 1) 

2 

when two indices are equal (type II), and the two sets of three states (m — —1, 0, 1) 



^ E ('^123)' ^nns , (4.9) 



*iT2f«3 ^= 71 E 6-^'="-^/^ (ai23)' (1 ± '^12) %nn.n, , (4.10) 



fe=0 



if all (ni, 77.2, 71.3) are different (type III). In this new orthonormal set of irreducible basis states 
the Fock representation of the Hamiltonian Hq reads 



The basis states \E'^''" are eigenfunctions of 0-123 and cri2T 

ai2TM'Jr'^ = ivff,™^^ . (4.11) 
Under ai2 and T separately, however, they transform into each other 

N N 
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We therefore have the following irreducible representations. The singlet states ^'('')+, the "axial" sin- 
glet states the doublets and the "axial" doublets («'(+^)-; ^'^-i)"). Since the 
partner states of the doublets transform into each other under the symmetry operations cri2 or T, the 
corresponding values of the energy functional are equal. 

The energy matrix elements of the irreducible states can then be expressed in terms of the basic 
matrix elements as given in Appendix C. Due to this decomposition of the Fock space into the irreducible 
sectors, the variational approach allows us to give upper bounds for states in each sector. The values 
of the energy functional for the states in each irreducible sector with the smallest number of knots 
= 4.8676 .g2/3, f [^^qo^^] = 7.1915 g^'^, ^[^012^] = 13.8817 g^/a^ and f [^^^2^^] = 15.6845 y^/s 
give first upper bounds for the lowest energy eigenvalues of the singlet, the doublet, the axial singlet, and 
the axial doublet states. 

In order to improve the upper bounds for each irreducible sector, we truncate the Fock space at a 
certain number of knots of the wave functions and search for the corresponding states in the truncated 
space with the lowest value of the energy functional. We achieve this by diagonalizing the corresponding 
truncated Hamiltonian fltrunk to find its eigenvalues and eigenstates. Due to the orthogonality of the 



truncated space to the remaining part of Fock space the value of the energy functional (4.2) for the 
eigenvectors of -fftrunk coincides with the i?trunk eigenvalues. 

Including all states in the singlet sector with up to 5 knots we find rapid convergence to the following 
energy expectation values for the three lowest states ^i, 6*2, S'a 

f [S'l] = 4.8067 g^l^ (4.8070 5^/^), 
E\S-^ = 8.2515 g^l^ (8.2639 5^/3), 

E\S-i\ = 9.5735 g^l^ (9.6298 5^/3)^ (4.12) 

where the numbers in brackets show the corresponding result when including only states up to 4 knots 
into the variational calculation. The lowest state £"1 , given explicitly as 

Sx ^ 0.9946 + 0.0253 - 0.0217 '^n'^J^ ~ 0.0970 "^nf-^^^ 

-0.0005 -^'-qS^ ~ 0.0033 ^-1,1^2^ - 0.0146 ^-^"^ - 0.0005 *^04^ 
+0.0040 «'[,°)5+ - 0.0080 ^-^20^ - 0.0038 ^-^1^2^ + 0.0001 1'[,°]5+ 

-0.0004 *[,"^4+ + 0.0011 ^-Jjalj^ - 0.0004 ^/f^^^ + 0.0031 *^2i^ , (4.13) 

nearly coincides with the state ^000^ > the contributions of the other states are quite small. Similarly 
including all states in the doublet sector with up to 6(5) knots the following energy expectation values 

for the three lowest states d[^^\ D2^^\ D^^^^ 



^[Dpi)] = 7.1682 5^/3 (7.1689 g^^^), 
^[^r^'] 9-6171 5^/^ (9.6394 .g^/^), 

£[D^^^^] = 10.9903 3^/3 (10.9951 5^/3). (4.I4) 

have been obtained. Including all states in the axial singlet sector with up to 8(7) knots we find the 
following energy expectation values for the three lowest states Ai,A2, A3 

£[Ai] = 13.2235 g^^^ (13.2275 g^^^), 
£[A2] = 16.6652 g^/s (16.7333 g^^^), 

SlAs] = 19.1470 g^/^ (19.3028 g^^^). (4.15) 

Finally taking into account all states in the axial doublet sector with up to 8(7) knots we find the following 
energy expectation values for the three lowest states c[^^\ C2^^\ 

Eld"^^^] = 14.8768 5^/3 (14.8796 g^^^), 
£[C^±i)] = 17.6648 5^/3 (17.6839 g^^^), 

f [C^±i)] = 19.9019 5^/3 (19.9914 g^/^)^ (4.I6) 
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We therefore obtain rather good estimates for the energies of the lowest states in the spin-0 sector. 
Extending to higher and higher numbers of knots in each sector we should be able to obtain the low 
energy spectrum in the spin-zero sector to arbitrarily high numerical accuracy. 

In summary comparing our results for the first few states in all sectors, we find that the lowest state 
appears in the singlet sector with energy 

£;n = 4.8067 5^/3 ^ (4.I7) 



with expected accuracy up to three digits after the dot. Its explicit form is given in ( 4.13| ) to the 



accuracy considered. For comparison with other work we remark that due to our boundary condition 



(3.9) all our spin-0 states correspond to the 0^ sector in the work of |16|] where a different gauge invariant 
representation of Yang-Mills mechanics has been used. Their state of lowest energy in this sector is 
9.52 g^/"^. Furthermore in using an analogy of SU{N) Yang-Mills quantum mechanics in the large 
A'' limit to membrane theory, obtain the energy values 6.4690 g^^^ and 19.8253 g^^^ for the groundstate 
and the first excited state. 



The expectation values for the squares of the electric and the magnetic fields for the groundstate ( [4.13D 
are found to be 

{0\E'^\0) = 6.4234 g^/s^ {0\B'^\0) = 3.1900 g^/^ , (4.18) 

and the value for the "gluon condensate" is therefore 

(OIG^IO) := 2 {{0\B^\Q) - (0|£;^|0)) = -6.4669 5^/3 . (4.I9) 

These results are expected to be accurate up to three digits after the dot. Hence the variational calculation 
shows that the vacuum is not self(anti-)dual and that a nonperturbative "gluon condensate" appears. 



B. Higher spin states 



For the discussion of the eigenstates of the Hamiltonian Hg^o with arbitrary spin we write 

Hg^O = Ho + Hspin (4.20) 



with the spin-0 Hamiltonian (4.1) discussed in the last subsection and the spin dependent part 



1 ^ 



1=1 



(a;2 - xlY 



i,j,k cyclic 



(4.21) 



Introducing t he lo wering and raising operators ^± ^1 ± 1^2 , the spin dependent part i/spin of the 
Hamiltonian (4.20) can be written in the form 



spm = 8 (^l + ) (Vl -V2) + - (e+C- + e-e+) (^1 + V2) + -esVs 



(4.22) 



Since the Hamiltonian (4.20) commutes with and J^, the energy eigenfunctions can be charac- 

terized by the two quantum numbers J and M. Furthermore we shall expand the wave function jm in 
the basis of the well-known D functions [Q, which are the common eigenstates of the operators — 
Jz and ^3 with the eigenvalues J, M and k respectively, 



I '2, J I 1 

*jM(a;i,a;2,a;3;a,/3,7) = ^ J —^--^^> jMk{xi, X2, X3)Dl:2j{a, P,"f) , (4.23) 



k=-.J 



where {a,(3,'y) are the Euler angles. We have the relations 



kM 



M ■ 



(4.24) 
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The task to find the spectrum of the Hamiltonian ( 4.2C| ) then reduces to the following eigenvalue problem 
for the expansion coefficients ^'jma: for fixed values of J and M 



E 

k=-J 



JMk 



= 



(4.25) 



Since the spin part i?spin of the Hamiltonian does not commute with ^3, nondiagonal terms arise, 
coupling different values of k. We shall in the following limit ourselves to the case of spin-1. Using the 
linear combinations 1451 



vl'i(a;i,a;2,a;3) := --j= ['i' j=i,M,k=i{xi, X2, X3) - '^j=i,M,k=-i{xi,X2,X3) 
^'2(a;i,a;2,a;3) := [^'j=i,M,fc=i(a;i, 2:2, 2^3) + '^j=i,M,k=-i{xi,X2,X3) 

^3{xi,X2,X3) -.^ ^J=l,M,k=o{xi,X2,X3) , 



(4.26) 

(4.27) 
(4.28) 



the corresponding eigenvalue problem ( 4.25 ) for spin-1 decouples to the following three Schrodinger 
equations for the wave functions '^aixi, X2, X3) 



2^ dx: 

2 — 1 



^aix) ^ E^aix) , a =1,2, 3 



(4.29) 



with the effective potential 



1 



1 



V:"{xuX2,X3) :=-(H + K) = 7T 



xl + xl 



xl +xl 



2\{xl-xir {xl-xlY 



a, b, c cyclic 



(4.30) 



In the spin-1 sector we have therefore succeeded to reduce the Schrodinger equation to three effective 
Schrodinger equations for the scalar degrees of freedom with an additional effective potential induced by 
the rotational degrees of freedom. Since the effective potentials V^^ are related via cyclic permutation 



C^123V^ =¥2 0-123 , 0-123^2° — ^3 '^123 , Cri23V3'' = V]° (7123 



(4.31) 



all energy levels in the spin-1 sector are threefold degenerate. 

As in the spin-0 sector we may use the variational approach to obtain an upper bound for the lowest 
spin-1 state. The variational ansatz 



*a(a;i,X2,X3) := (.t^ - xg)(x^ - .t^) J| ^'o(a;„ a;,) 



(4.32) 



satisfies both the boundary conditions (3^) and ( 3.10 ) and vanishes at the singularities of the additional 
effective spin-1 potential Vcg. For the optimal values 



LOa = 1-1814 5^/3, uJb^LUc^ 2.34945 g^^^ , 
we obtain the energy minimum 

-B,pin-i = 8.6044 5^/3 . 



(4.33) 



(4.34) 



Analogous treatments of higher spin states can be carried out correspondingly. Using the linear combi- 
nations [E5i 



1 



^j^l^^{xi,X2,X3) —j=['^J^M.k{^1^^2,X3) ±'^J^AI^^k{xi,X2,X3)] , fc^O 

*jo(a;i,a;2,a:3) := *j,j\/,fc=o(a^i, a;2, 2:3) , 



(4.35) 
(4.36) 
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and noting that there are no transitions between the states ^'j^-i with even and odd k, and with + and 
— index, the corresponding eigenvalue problem (4.25) for spin-J decouples into four separate Schrodinger 
eigenvalue problems. For spin-2 one finds one cyclic triplet of degenerate eigensates and two singlets 
under cyclic permutation, for spin-3 two cyclic triplets each consisting of three degenerate states an d one 
singlet, and so on. The corresponding reduction on the classical level using the integrals of motion (3.11) 

has been done in |11|. 

We conclude this subsection by pointing out that our variational result (4.34) shows that the higher spin 
states appear already at rather low energies and therefore have to be taken into account in calculations 
of the low energy spectrum of Yang-Mills theories. 



V. CALCULATION OF THE TOPOLOGICAL SUSCEPTIBILITY 



The explicit evaluation of the Witten formula (3.16) for the topological susceptibility allows us to check 



the consistency of the results for the low energy spectrum obtained in Section IV using the variational 
approach. 

Using the groundstate Si in ( 4.13| ), obtained from minimization of the energy functional in the singlet 



sector including irreducible states with up to 5 knots, and the expressions for the matrix elements of 
in the basis of irreducible states given in Appendix C, we obtain 



contact 



(01 f S2|0) = +0.0005117 .g"/3 (+0.0005119 g^^/^) (5.1) 
v zvr/ 



for the contact term in the Witten formula. The number in brackets gives the corresponding result for 
up to 4 knots. 

Since the Q-operator is a spin-0 operator and symmetric under cyclic permutations, the propagator 
term involves only the singlet states in the spin-0 sector. Using the formula for the matrix elements of the 
topological operator Q stated in Appendix C and including the lowest fifteen (ten) excitations S2, ■ ■ ■ , Siq 
(S2, . . . , 5*11) obtained approximately in the variational calculation as eigenvectors of the truncated Fock 
space including irreducible singlet states up to 5 knots (4 knots), we obtain |^ 

^-^in = -2y = -0.0004819 g''^' (-0.0004622 g'^^^) . (5.2) 



We see that the sum of the contact contribution (5.1) and the propagator contribution (|5.2|) seem to 



tend to zero when extending the variational calculation to Fock states of higher and higher number of 
knots. For comparison w e po int out that using the irreducible sing let states ^'ooo^ , ^'ooi"^ , ... up to 5 



knots (4 knots) in (4.8)-( 4.10 ) directly, instead of the eigenstates Si, S2, ■ ■ ■ , Sie (Si, S2, ■ ■ ■ , Sn), we 
get -0.0005205 ^"/^or the contact contribution (|1]) and -0.0003808 g^^/^ (-0.0003761 g^^^^) for 
the propagator contribution ( ^.2| ). We herewith find stro ng support that our variational results are in 
accordance with vanishing topological susceptibility ( [3.16[ ). 



VI. CONCLUDING REMARKS 



In this paper we have analysed the quantum mechanics of spatially homogeneous gauge invariant SU{2) 
gluon fields with theta angle. We have reduced the eigenvalue problem of the Hamiltonian of this toy 
model for arbitrary theta angle to the corresponding problem with zero theta angle. The groundstate, 
its energy eigenvalue, its magnetic and electric properties, as well the corresponding value of the "gluon 



Here the lowest six excitations S2, ■ ■ ■ ,87 are found to give the contributions, -103.3 • 10"" g"/^ (-107.7 • 
lQ-(i giifi)^ -201.6 • IQ-^ <?"/^ (-205.3 • 10"'' g"''^), -124.1 ■ 10"" <7"''^(-120.4 ■ 10"" ff"/^), -8.8 ■ 
10-6 <;i4/3(_9 3 . 10-6 gi4/3)^ _27.3 ■ 10"" ffi4/3(-18.4 • 10"" ff"''^) _o^g . ^g-e ^14/3 . ^g-e ^i4/3) 

respectively. The contributions from the remaining higher excitations Ss, ■ ■ ■ , Sie (Ss, ■ ■ ■ , Sn for up to 4 knots) 
are of the order of 5 ■ 10"" or less and form a series which is rapidly decreasing with the number of knots. 
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condensate" and the lowest excitations have been obtained with high accuracy using the variational 
approach. Furthermore it has been shown that higher spin states become already relevant at rather low 
energy. The groundstate energy has been found to be independent of the theta angle by construction 
of the explicit transformation relating the Hamiltonians with different theta parameter. An explicit 
calculation of the Witten formula for the topological susceptibility using our variational results for the 
groundstate and the low lying excitations give strong support for the independence of the groundstate 
energy of theta thus indicating the consistence of our results. We have found a continuous spectrum 
and the corresponding eigenstates of the topological operator in this approximation and shown that its 
characteristics coincide with the Euclidean self(anti-)dual zero-energy solutions of the classical equations 
of motion. They are the analogs of the instanton solutions, but do not correspond to quantum tunneling 
between different vacua. The generalization of these investigations to SU (2) field theory following is 
under present investigation. 
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APPENDIX A: TOPOLOGICAL CHARGE OPERATOR, ZERO-ENERGY SOLUTIONS OF 
THE CLASSICAL EUCLIDEAN EQUATIONS OF MOTION AND TUNNELING 

AMPLITUDES 

In this Appendix the solution of the eigenvalue problem for the topological charge operator is described 
and the relati on be tween its characteristics and the Euclidean zero energy trajectories of the unconstrained 
Hamiltonian ( 2.25| ) discussed. We shall also discuss the role of these Euclidean zero energy solutions of 



the classical equations of motion to tunneling from one valley to another. 



1. The eigenvalue problem for the Q- operator 



The eigenvalue problem for the topological charge operator 

g|vi/(t))^ = A|vi/(t))^ (Ai) 

in the Schrodinger representation reduces to the solution for the following linear partial differential equa- 
tion 



XiX2-^^\{xi,X2,X3) + X2X-i-^'^\(xi,X2,X3) + X'iXi-^'^ \{xi , X2, X^) = -i A* A (xi , X2 , X3) . (A2) 

0x3, 0x1 0x2 



The conventional method of characteristics relates this problem to the solution of the set of ordinary 
differential equations 

dxi dx2 dx3 ^^^^ 



X2X3 X3X1 X2X1 



The integral curves corresponding to (A3) can be written in the form 

h = x\ -x\ , h^xl-xl . (A4) 
These integral curves promt us with the introduction of the new adapted coordinates (C, f], p) 

C := a;i , 'q:=xl-xl, p:=xl-xl. (A5) 
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Such functions can be used as a suitable coordinates on the subset 

< < 2:2 < < cx) 



(A6) 



of the whole configuration space . The subset (A6) corresponds to the domain < C < Vv + C^ < 
+ C^- -Due to the symmetry of the Q -operator under arbitrary permutations of the canonical pairs 
Xi,pi the results can be extended to the whole . 
Writing the wave function in terms of new variables 



the partial differential equation (|A2|) reduces to the following ordinary differential equation 



d 



Ve + Wc^ + Pg^wx (c, V, p) = ~i\^wx{(:, V, p) 

The general solution of this equation can be written in the form 



W\{C,,'n,p) = 1'o(?/,p) exp[iA-^^F farctan 

9 Vp \ \y/n 



p-v 



(A7) 



(A8) 



(A9) 



with the arbitray function 5*0 (?7, p) and the Jacobi elliptic integrals F{z, k) of the first kind |4^. In terms 
of the original coordinates {xi,X2,X3) the eigenfuctions for the topological charge operator in the sector 
Xi < X2 < x^ therefore have the form 



^\{xi,X2,X3) ^ 'i'oixl - xl,xl - xl)exp 



iX- 



87:' 



zF arctan 



Xl 



V4 











X2 






— x\ ^ 


) 



(AlO) 



In the other sectors the corre sponding wavefunction is obtained from ( AlC ) by cyclic permutation. Note 
that the eigenfunctions ( AlO ), which constitute the most general solu tion o f the eigenvalue problem (Al) 
for the Q operator, do not satisfy the boundary conditions ( [3.S| ) and ( 3.10| ) necessary for the Hermiticity 
of the Q operator. 

In the next paragraph we will show that the characteristics of the topological charge operator ( |A3| ) coin- 
cide with of the equations which determine the zero-energy solutions of the Euclidean classical equations 
of motion. 



2. Euclidean zero energy trajectories in Yang-Mills mechanics 

The Euclidean action 

is obtained from the corresponding action in Minkowski space by inverting the potential V{x) 
In the one dimensional case the solutions of equation 

dx 
d7 



coorespond to trajectories with zero Euclidean energy 

1 



= -x' ~ V{x) , 

and the same time satisfies the classical Euclidean equations of motion 

'x 

dr'^ dx 



(All) 
-V{x) . 

(A12) 
(A13) 
(A14) 
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Such type of trajectories play an important role in the description of quantum mechanical tunneling 
phenomena [^7| . In the case that the potential V{x) has at least two local minima, say at a; = —a and 
X — +a, with V ^ 0, the Euclidean zero energy trajectories starting at —a and ending at a correspond 
to quantum tunneling into the classically forbidden region. The Euclidean action for these classical 
^Euci ^ trajectories 



:rEucl I 



l_E=0 



= / dt 



dxy/2V{x) , 



(A15) 



determine in the semiclassical limit the WKB amplitude for a particle to tunnel from x = —a to x — +a 
\T{E = 0)1 = cxp 



dxy/2V{x) 



{l + 0{h)) 



(A16) 



The potential of the unconstrained system considered in the present article has three valleys. The 
question is whether there exist the trajectories corresponding to tunneling between the valleys. To 
answer the question let us rescale the coordinates Xi — > g~^Xi and write down the Euclidean action of 
the model in the form 



^Eucl 



1 

2? 



dr ^ (i^ + x^x^ 



cycl. 



The equations of motion then read 

Xi — Xi{x'^ + x1); J, j, A: cyclic . 
The class of trajectories with zero energy 



•2 2 2 

X,X 



cycl. 



can be chosen as the solutions of the following system of first order equations 

±1 = ±X2X3 , ±2 = ±X3Xi , ±3 = ±X\X2 ■ 



(A17) 



(A18) 



(A19) 



(A20) 



If we fix one and the same s ign on the r.h.s. of all equations ( A20 ) the n th ey completely coincide with 
the characteristic equations ( [A3| ) of the Q-operator. Furthermore from (A4) we see that the zero energy 
Euclidean solutions admit no trajectories from one V = minimum to another and they have no relation 
to the quantum tunneling phenomena. 



3. Q- operator and selfdual states 



The commutator of the Hamiltonian of (4.1) in the spin-0 sector and the topological charge Q 



9 

[Ho, Q] = {{P1P2X3 + P2P3X1 + P3P1X2) - gxiX2X3{xl +xl+ xj)) (A21) 

vanishes only in the subspace of states j^* > which satisfy the Euclidean self(anti-)duality conditions 

PjI^-) = ±gXjXk\'i'); k cyclic , (A22) 



which are the quantum analogs of the Euclidean E = constraints (A20) discussed before. Rewriting 
the Hamiltonian Ho in ([4. if) in the form 



1 Stt 
Ho = - iPi+ 9^^^j^k) = ^P'^ gXjXkf ± — ( 

9 



(A23) 



i.j,k cycl. 



i.j,k cycl. 



we see that the Hamil tonia n Ho and the topological operator Q coincide on the subspace of the Euclidean 
self(anti-)dual states (A22) 
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^I*> = T%Qm ■ (A24) 
9 

Comparing this discussion with the corresponding original situation H = (1/2) (H"^ + Bf^) and Q — 
~{as/2TT)Il°;B°' in terms of the constrained fields 11° and Af, where Hq — / 9'^)Q only in the 

Euclidean self(-anti)dual case Ilf = Bf, we see that Eq. ( [A22| ) correspond to the unconstrained analogs of 
the self(anti-)dual configurations in Euclidean space. The question arises whether there are any self(anti- 
)dual states (A22) which are both eigenstates of Q and H. The solution 

*gj3 := exp[^igxiX2X3] , (A25) 

of the self(anti-)duality conditions ( A22| ) in Euclidean space is neither eigenfunction of Hq nor of Q 

^o*SD - T^g^to = ±2^*fD > (A26) 



The well-known exact nonnormalizablc zero-energy solution of the spin-0 Schrodinger equation (4.1), 

^ Aexp[TgxiX2X3] , (A27) 

which differs from functional up to a factor of i in the exponent, actually satisfies the self(-anti)duality 
conditions in Minkowski space 

■A 

dx_ 

(corresponding to = ±ii?f ), but is not an eigenfunction of the topological operator Q 
Stt 
9 



i—-^t = ^igxjXk-^t; hj, k cyclic , (A28) 



= ±ig^{xlxl + xjxl + xjxD-iiQ = ±2iV{xi, X2, X3)'^q ■ (A29) 



Finally we point out that the (approximate) groundstate wave function (4.13) obtained in the variational 
approach is not self(anti-)dual. 



APPENDIX B: LOWER BOUND FOR THE SPIN-0 HAMILTONIAN Ho 



In this appendix we would like to derive a lower bound for the spin-0 Hamiltonian Hq in (4.1) along the 



line of Ref . |lj] . Using the boundary conditions 
inequality for oscilator on positive half line 

dx'^ 

it follows that 



|3.9| ) and ( p. 10 ) and based on the well-known operator 



Ho>2 (-^ + sV2g{xi +X2+ xs)) =: ^H' 



(Bl) 



Since the Hamiltonian H' is known to have a discret spectrum, this is true also for Hq. An important 
open question is at which energy the groundstate is. The knowledge of the groundstate energy of H' 
in inequality (Bl) would provide a lower bound for the groundstate energy of Hq. Due to the additive 
structure of the potential term in H' one can make a separable ansatz for the solution of the corresponding 
eigenvalue prob lem. The e nergy of the lowest such separable H' eigenstate satisfying the above boundary 
conditions (3^) and ( 3.10 ) is 



F' 

scp 



6|eo|(35/2)'/'- 9.1924 g^/s ^ 



(B2) 



where ~ —2.3381 is the first zero of the Airy function. From the operator inequality (Bl) and the lower 
bound ( ^^ ) for separable solutions of H' we therefore obtain the lower bound of the energy functional 
for separable functions 



1 



f [*scp] > TT^^scp = 4.5962 5^/3 



(B3) 



Finally we remark, as has been pointed out already in [ pl| , that an analogous variational calculation for 
H' sho ws th at also the groundstate energy of the Hamiltonian H' in Eq. (Bl) is lower than the value 
iJgpp of (B2) for the lowest separable solution. 
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APPENDIX C: MATRIX ELEMENTS 



For the evaluation of the energy functional, the calculation of the value of the " gluon condensate" of 
the groundstate, as well as the propagator term in the Witten formula we need the matrix elements of 
E^^B^ and Q with respect to the irreducible Fock space states ( |4.8| )-(4.1C) built from the basic Fock space 
states (O). 



1. Basic matrix elements for Hamiltonian and topological charge 



The matrix elements of E'^ and B'^ with respect to the basic Fock space states (4.7) with (cji = uj2 
CJ3 — 3^/"^ 9^^^) are given by 



(^mi;iTi2;m3 1^ I^ni;n2;n3) ~ ^ 9 ^ ^ ^ '^mini^rrijrijSn 

cyclic 



and 



where 



mi ;m2 ;m3 



cyclic 



^mn (5m«(2rt + 3/2) ± 5m{7i+i) \/n + 3/2Vn + 1 ± (5(m+i)„\/n+ 1/2^" ■ 
For the topological operator Q we have 

2^58/3 



IQI 



^7/231/6 



V Q" Q+ Q+ 

1/6 / y ^miUi ^mjiij ^irikTik 



cyclic 



where 



+ _ 1 (-l)"+"(2m+ l)!!(2n+ 1)!! 

l-4(m-n)2 ^(2m + l)!(2n^!y! 



m-n (-l)"+"(2m+ l)!!(2n + 1)!! 



1 - 4(m - n)2 ^(2m + l)!(2n+ 1)! 



(CI) 
(C2) 

(C3) 
(C4) 

(C5) 
(C6) 



2. The irreducible matrix elements in terms of the basic ones 



For any operator O invariant under the permutations ct^ , such as E'^, i?^, the Hamiltonian Hq and 
the topological operator Q, 



(C7) 



t he m atrix elements of the irreducible states (5'^^''^|0|5'^''^) of type I (4^), type II ( |4.9| ), and type III 
(4.10) can then be expressed in terms of the basic matrix elements 



■' 'nin2n3 \ * irii ;m2 ;ni3 I "-^ I * "l ;n2 ;n3 , 

as follows. For the type I, II and III singlet states we have 



jmmm 



(C8) 
(C9) 



(CIO) 



and 
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'n2n3"l 



respectively. The transition elements between the type I, II and III singlets are 

(*(^);L|0|vi/(°l+) = V3A^™ (C12) 
(*(^H,|0|*(°)+„3) = x/6A^~3 (C13) 

and 

(^(0)+ |C)|^r(o)+ ) = 727^™""^ + \/27W™™'' + \/2A^'"""' . (C14) 
For the doublets, which exist only for type II and III states, we have 

/^(1.2)+|Q|^(l,2) + \ ^ ^7n7nr _ ^7n7nr /(^-^qn 
\ mmr 1^1 nns I ^nns ^nsn \^ ^ / 



and 



|(-)|iTfCl,2)+ \ _ A^mimsma _ /i /o\ A^mim2m3 _ /i A^mim2m3 

\ mi 7712 7713 l l ni 712 713/ *■ 71 1 712 713 V / /*^ '-7l371i7l2 \ / /*^ '-7l2 7l371i 

+-M"2Cnr - (i/2)A^:r3\"r' - (i/2)A^:rM:r^ (cie) 



for the type III doublets. Their transition elements are 

/^(l,2)+|Q|^(l,2)+ \ /2^77i77ir _ |'^7n77ir ^ ^7n77i,r )/^ (QIJ) 

\^ mmr 1^1^711712713/ ^ '^■'"^'77.1712713 V n3"in2 '712713711 y / ^ ^ ■ \^ ^ ' J 

For the axial singlets we have 



^71117112 7113 
'7I27I37I1 



(\t/(0)- \0\\I/(°'>- ) = /M'niiTiama , A|mi7ii2m3 1 yW' 

\ 711171127113 I ^ I 7li7l27l3 / '71i7l27l3 ' '7I37I17I2 ' 7 

A^71li71l271l3 A y< 711 1711 2 711 3 A "l 1 7112 71l3 ^'P1S^ 

-■^^772771773 " "'^773772711 " -^^77 1 713 772 ' {'^ ) 

For the axial doublets we have 

/\I>(1'2)- 1/01^(1,2)- \ _ yW717im27l73 _ ( -I /0^ yW "1 1 7112 7173 _ f 1 /0^ yW711l71l2m3 
\ 77li71l2 77l3 I ^ I 77 1 77 2 713 / '77i7l2773 V / Z*^ '7l371i772 V / Z*^ '7I27I37I1 

-■M"2C«r + (i/2)-M"3nX"^ + (i/2)A<™r^ • (C19) 
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